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Abstract 

We consider a 0-ir Josephson junction consisting of asymmetric and tt regions of different 
lengths Lq and L n having different critical current densities j c o and j Ct7T . If both segments are 
rather short, the whole junction can be described by an effective current-phase relation for the 
spatially averaged phase tp, which includes the usual term oc sin(^), a negative second harmonic 
term oc s'm(2ip) as well as the unusual term oc Hcosip tunable by magnetic field H. Thus one 
obtains an electronically tunable current-phase relation. At H = this corresponds to the <p 
Josephson junction. 

PACS numbers: 74.50.+r, 85.25. Cp 
Keywords: ip Josephson junction 



C/3 



1 



I. INTRODUCTION 



Recently we proposed^ to implement a tp Josephson junction (JJ)2 with magnetic-field 
tunable current-phase relation (CPR) based on an 0-ir JJ with the and 7r segments of differ- 
ent length Lq 7^ L n . This proposal was made keeping in mind YB^CuaCVNb ramp zigzag 
JJ technology^ (or a similar one^ with Nd2- x Ce x Cu04-Nb) established recently in our 
group also^. However, in experiment we were more successful^ in employing superconductor- 
insulator-ferromagnet-superconductor (SIFS) O-ir JJs^~— , where the lengths of and n seg- 
ments are equal, but critical current densities j C) o and j C 7r in the and 7r parts are different. 

Therefore, in this paper we present a more general theory, which describes an effective 
tp JJ made of asymmetric and n regions of different lengths L and L n having different 
critical current densities j Cj o and j Ci7r . 

II. MODEL 

The static sine-Gordon equation that describes the behavior of the Josephson phase (ft in 
a 0-7T JJ is 



Here /io is the magnetic flux quantum, n^dj is the specific inductance (per square) of the 
superconducting electrodes forming the JJ and j is the bias current density. The prime 
denotes the partial derivatives with respect to coordinate x. We assume that the critical 
current density j c (x) has the form of a step-function 




-j c (x) sin0 = -j. 



(1) 



jc = jc,0 > 



< x < U 



-Ltt < X < 0. 



(2) 
(3) 



We write the critical current density j c (x) as 



jc(x) = (jc(x)) \l+g(x)}, 



(4) 



where 
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is the average critical current density, L = Lq + L n is the total length of the junction, and 
(g(x)) = 0. The function g(x) is defined as 

9W = Im> - 1 (6) 

that results in 

, . l9o, < x < Lq, 
9{x) = I (7) 

yg w , —L n < x < 0. 

where 

(jc,0 — jc,ir)Ln (jc,0 — jc,n)L , . 

fi 10 = ~ — F — r~ — r - ' ^ ~~ — : — f — r~ — r~ ■ W 

,7r-^7r 7T 

Then we divide Eq. (JTJ) by | (j c ) | and normalize the coordinate x to the Josephson length 
calculated using | (j c ) |, i.e., 

Xj = V 2^] (J C ) | • (9) 

Thus, we obtain a normalized sine-Gordon equation for the phase difference (f>(x) 

<p" - sgn(Oc)) [1 + fi 1 ^)] sin0 = -7, (10) 

where 7 = j/| (j c ) I is the normalized bias current density. It is worth mentioning that (j c ) 
can be positive as well as negative. Below, for the same of simplicity, we assume (j c ) > 0. 
Thus, Eq. ffTQ]) becomes 

<p" - [l+#(x)]sin0= -7. (11) 

In the case (j c ) < the substitution — > tt — converts Eq. ( flOi) to the same Eq. ( 1TTT) . 
We look for a solution of Eq. ( HID in the form 



0(x) = -0 + £(x) sin-0, (12) 

where 

^ = m*)) (13) 

is a constant average phase, while £(x)sin-0 describes the deviation of the phase from 
the average value, i.e., (£(x)) = 0. Further we assume that the deviation is small, i.e., 
|£(a;) sin^l <C 1. Then we plug the relation ffl2l into Eq. (fTTj) . expand it in series in 
sinijj, and keep the terms of zero and first order. We get 

siiaip — [1 + gf(x)][l + cosip] sinip = —7. (14) 



The constant terms (zero order of £ in Eq. (1141) ) are 

7 = sixit/j + (£(x)g(x)) cos ^ sin -0. (15) 

The terms of first order of in Eq. f[T4l) are 

T - <7(z) = {£ + - (t(x)g(x))} cos^. (16) 

Numerical calculations show that the two terms oc cosip have an extremely weak effect 
on solutions of Eq. (TIB"]) . We neglect these terms and obtain for £(x) 

e - g(x) = o. (17) 

We treat solutions of Eq. (|T7|) by using the matching continuity (at x — 0) and boundary 
(at x = —l T = —L n /Xj, x — Iq = L /\j) conditions 

Gr(o) = G>(o), &(o) = £(o), (is) 

sin?/; = h, &(Z ) sintp = h. (19) 
The applied field H is normalized by H c \/2, i.e., 

h = 7T- H « = JST- < 20) 



where A is the effective magnetic thickness of the JJ. We integrate Eq. ( fTTj) once and obtain 

&(x) = g (x - lo) + -r^-r, < x < / , (21) 
smip 

C(x) = g n (x + Q + < x < 0. (22) 

sm-0 

The second integration results in 

= flb f T " /oX ) + ^7 + ^ ( 23 ) 



for < x < / 



0i 

&r(aO = ^ f %r + ^ + -^7 + C, (24) 
for — l w < x < 0. 

The integration constant C can be obtained using the condition (£(#)) = 

C = ^{ 3 irn^J- (25) 



We use Eqs. (|23|) . and ( 124)) and obtain the average (£,(x)g(x)) in the form 

(£(x)g{x)) = r + r h -\, (26) 
where the coefficients T and are given by 

72/2 (a _ t ^2 
-p t < '7r Uc,0 Jc,7rJ Cy7\ 

0_ O 7""* 1 T~ T^2' ^ 

3 U c ,oto + 3c, Jn) 

■p, Jc,0 Jc,7T /r)0\ 

ih = ^ — : — 7 — r~ — r- ^> 

Using Eqs. ( 1T51) and ( 1261) we find the current-phase relation in the form 

j — (jc) (sin if) + T sin ip cos -0 + hY h cos -0) . (29) 

It is worth noting that there is a simple relation between the coefficients To and 1~V 
Indeed, it follows from Eqs. (ETJ) and Q2B} that 

r = -^r 2 . (30) 

In the case of equal lengths of amd tt parts (Iq = l w = 1/2) we find 

r = -— ( i'>o-i'A \ r h = . (3i) 

The energy [/ (?/>) corresponding to the current-phase relation (1291) is given by 

U{$) = (j c ) (l - cos ^ + M\ sin ^ + y sin 2 J . (32) 



III. CONCLUSIONS 

We have extended our previous results^ to the case of arbitrary critical current densities 
jc,o 7^ jc,n more relevant for experiment^. The dependence (1291) of the CPR on the phase 
and applied field is the same as in our previous studyi. The difference is in the formulas 
$ZE]) for T and r h . 
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